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Abstract 

We describe how a generalization of the original Gross-Neveu model from 
U{N) to 0{N) flavor symmetry leads to the appearance of a pairing conden- 
sate at high density, in agreement with the conjectured phenoinenon of color 
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superconductivity in (3 + l)-dimensional QCD. Moreover, the model displays 
a rich phase structure which closely resembles the one expected in two-flavor 
QCD. 



I. INTRODUCTION 

The work that will be described in this note had its origin at the beginning of the 
history of TMU-Yale collaboration about a decade ago, when one of us (A.C) was working 
on QED with a chemical potential ||1],0] and another (H.M.) suggested that the study of 
the Gross-Neveu model with a chemical potential might have interesting applications 
to doped polyacetylene. This resulted in a pair of papers, in the first of which 0] the 
effective potential for the Gross-Neveu model was computed with a non-vanishing chemical 
potential fi at leading order in large N. In agreement with earlier work p, a first-order 
phase transition was found at /i = m/-\/2, where m is the fermion mass, and this result was 
successfully compared to the phenomenology of doped polyacetylene [H,§. In the second 
paper [0], the thermodynamic Bethe ansatz was employed to compute the next-to- leading 
1/A^ correction to the position of the phase-transition. 

But the Gross-Neveu model, possessing aymptotic freedom and spontaneous chiral sym- 
metry breaking, was originally introduced as an analogue of QCD. With renewed interest 
in QCD at finite density |^, the TMU-Yale collaboration shifted its focus to trying to un- 
derstand the phenomenon of color superconductivity. Before describing that work, we make 
a few remarks about the role played by the large N approximation, which will be used 
throughout this paper, in the prediction of the superconducting phase transition. 

The issue is the following: the superconducting phase will involve a non- vanishing pairing 
condensate, which spontaneously breaks the U{1) of fermion number. But the generalized 
Gross-Neveu model is a 1 -|- 1 dimensional theory, and theorems exist that forbid the 
spontaneous breaking of continuous symmetries in 1 -|- 1 dimensions. Nevertheless, to leading 
order in 1/A^ one can compute the effective potential, and one finds that the condensate 
indeed forms. As Witten has explained |T0[, corrections are expected beyond the leading 
order which will remove the condensate. More explicitly, the vacuum expectation value of 
B[x) must vanish at any finite A^ because its correlation function {B{x)B'^{y)) will fall off, 
for large (x — y), as | x — y |~^/^. If we take the limit A^ ^ oo first before taking the limit 
\ X — y \^ oo it implies non- vanishing vacuum expectation value of B{x). 



One can enhance one's intuition on this matter if one accepts the analogy between the 
infinite number of fields that are present in the large N approximation and the infinite 
set of fields one obtains upon making a Kaluza-Klein reduction from a higher number of 
dimensions. Admittedly, this analogy is not perfect, because in the large N approximation 
the fields are all identical, whereas in Kaluza-Klein they typically represent a tower of states 
with indefinitely increasing masses. Nevertheless, the analogy suggests that the leading N 
theory somehow shares the physics of a theory in a higher number of dimensions, where 
the no-go theorems do not apply. The analogy further suggests that if one is interested 
primarily in elucidating properties of (3 + l)-dimensional QCD, the large N approximation 
to the Gross-Neveu model may be a more appropriate object of study than the full theory. 

Even more relevant examples in this context are provided by the Eguchi-Kawai model 



Tl[ , or by a class of large-N matrix models [p!2|,p!3| . In these matrix models the eigenvalues of 



the matrix are somehow converted, in the large-N limit, into the spatial dimensions. Thus, 
the matrix models in zero dimension thereby serve as definitions of 4-dimensional Yang- 
Mills and higher-dimensional superstring or supergravity theories. It is not clear to us if 
the analogy goes through to the vector-like large A^ theories such as ours. Nonetheless, it 
is at least suggestive enough to lead to the speculation that Coleman's theorem is violated 
in (1 + l)-dimensional large N theories because they are actually the theories in higher 
dimensions. 

II. THE 0{N) GROSS-NEVEU MODEL 

The phase structure of the original Gross-Neveu model is well-known [^ : at low temper- 
ature and density, the chiral symmetry is broken; as either is increased, the chiral condensate 
vanishes and the symmetry is restored. But there is no evidence for a second, pairing con- 
densate that would be the analogue of color superconductivity in QCD. It is necessary to 
add to the original model an extra term in the Lagrangian which, as it turns out, has a 
simple interpretation. The original model, with A^ identical fermi fields, possesses U{N) 
symmetry. The extra term is what can be added if one demands only 0{N) symmetry in 



flavor space. 

Explicitly, our lagrange density reads 






The flavor indices i,j are summed from 1 to A^. The first 2 terms are the original GN model. 
The third term, because of the way the flavor indices are summed, is 0{N) but not U{N) 
invariant. The last term represents the chemical potential. The gamma matrices are 2x2, 
and we choose 7° = cri, 7"*^ = —ia2; 75 = 0-3. Then 

2G2^«75^^^'^ V'^Ss^^^'^ = -G' [e„;3^t{^ VJ^^^] [e,s^?4'^] ■ (2) 

Following standard techniques [|I4| we introduce auxiliary fields: 



m{x) = —g^'ipip 

B{x) = -Gh^pi^^^i,f (3) 

and integrate out the fermions to obtain the effective action Fejj as a function of m and B. 
This is described in more detail in refs. [0,|1^]. 



As appropriate to a vacuum solution, we assume that m and B are constants, and write 

Teff = -N j{Sx)V,ff. (4) 

We set g'^N = A, and G^N = k/4. The large N limit is defined by letting N —>■ 00 with A 
and K fixed. We find 

Veffim, M) = ^ + ^ + Vj}}{m, M) (5) 



where M^ = AB'^B. Explicit computation yields ||T6| 

v}jUm, M) = --^ r dk[k+ + k_ + hn{l + e-^''+) + |/n(l + e"^*^")] (6) 
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where A is an ultraviolet cutoff, (3 = 1/kT (T = temperature) and 



fc± = Jh ± 262 (7) 



with 



bi = M' + m' + fi' + k' (8) 



and 



h^ = [M^m^ + ^\k'' + m')Y'K (9) 

III. RENORMALIZATION 

As with the original GN model, we renormalize the theory, thereby absorbing the de- 
pendence on A in the relation between the bare couplings A and k and their renormalized 
counterparts A/? and k^r- Recall that in the original model, one trades the dependence on 
A for dependence on a renormalization scale ttiq, or equivalently on the fermion mass mp- 
Once one fixes mp at its physical value, there are no further parameters in the theory: the 
renormalized coupling Xr varies with ttlq in such a way as to keep nip fixed. 

Similar features appear in the more general model, except that now there are two renor- 
malized couplings. As a consequence, the theory has two physical parameters: the fermion 
mass mp, and another parameter 5 that is independent of the choice of renormalization 
scale, and is defined by 

- ^ ^ ^ ^ (10) 



kr 2\fi K 2A 



(Although K, and A separately depend on the cutoff A, the combination represented by 6 does 
not.) Further computation |T6| allows us to express the renormalized Veff as a function of 



the condensates m and M, the physical parameters rrip and 6, and the temperature T and 
chemical potential /i: 
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1 ^oo 2 

Veff = 6M^-— dk[k+ + A;„ + -(In (1 + 6"'^'=+) + In (1 + e^'^'^-)) - 2ki 
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/9 
A;2 + m% 



- 2k{m' + M') ^ =] . (11) 



This expression is valid for 5 > 0. For 5 < 0, there is no fermion mass, and instead 
at /i = T = there is a pairing condensate M^ = A^. Since QCD has broken chiral 
symmetry but no pairing condensate at /i = T = 0, we shall restrict ourselves to 5 > 
henceforth. By studying how A^j and /t^ vary with renormalization scale, one sees that both 
are asymptotically free. (Because 5 is independent of this scale, it follows that if one of these 
couphngs is asymptotically free so must the other be.) 

IV. PHASE STRUCTURE 

The phase structure of the theory can now be mapped out by fixing mp^ 5, /i and T, and 
finding the minimum of Veff with respect to m and M . If m = M = 0, the theory is fully 
symmetric; if m 7^ 0, chiral symmetry is spontaneously broken, while if M 7^ 0, there is a 
pairing condensate which represents the analogue of color superconductivity. 

Because of the complexity of Vg//, this analysis must be done numerically. The results 
are schematically displayed in Figure 1 for a typical value of 5. 
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FIG. 1. A schematic illustration of the phase diagram on T — n plane. 



The point x is a tricritical point, at which the order of the phase transition changes from 

2 

first (solid line) to second (dotted line). The point P is located at /i^ = -^[l — e^'^'^^]. This 
generalizes the earlier result for the U{N) GN model, //^ = ^m'jp, for which 6 -^ oo. Note 
that, interestingly, even though we have asymptotic freedom, the M ^ phase persists for 
arbitrarily large fi. In fact, the transition temperature from the Cooper pairing phase to the 
unbroken phase is independent of fj,. 
The figure is qualitatively valid for 

1.13 



S> Sr. 
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For smaller, positive 6, the tricritical point disappears and the phase transition for fixed n 
is always second order. For 6 < the chirally broken phase disappears, and there is only a 
pairing phase at low T and an unbroken phase at high T. 

There is no case of a phase in which both m and M are non-zero. However, there can be 
coexistence of the m ^ and M ^ phases along the line of first-order phase transitions. 



V. CONCLUSIONS 

The phase structure displayed in Figure 1 is strikingly similar to that conjectured in the 
literature for two-flavor QCD. This gives us confidence that further properties of QCD can 
perhaps be studied reliably in this 1 + 1- dimensional theory, whereas calculations involving 
the full 3 + 1-dimensional QCD may be prohibitively difficult. For example, work is in 



progress [Tj] on djTiamical calculations, in which the system is prepared initially in a dense, 
hot state and allowed to cool and expand. Various distributions, such as the "pion" wave 
function and the "pion" correlation function, are followed as this evolution takes place. So far 
these computations have been performed only in the U{N) model, and in the approximation 



that the expectation value of physical quantities depends only on proper time r = y/f^ — z'^, 
but more general results are expected in future investigations. 
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